A unifying model for chemical and colloidal gels 
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We investigate the slow dynamics in gelling systems by 
means of MonteCarlo simulations on the cubic lattice of a 
minimal statistical mechanics model. By opportunely vary- 
ing some model parameter we are able to describe a crossover 
from the chemical gelation behaviour to dynamics more typ- 
ical of colloidal systems. The results suggest a novel connec- 
tion linking classical gelation, as originally described by Flory, 
to more recent results on colloidal systems. PACS: 05.20.-y, 
82.70.Gg, 83.10.Nn 



The chemical gelation transition, as it is typically ob- 
served in polymer systems, transforms a solution of poly- 
meric molecules, the sol, from a viscous fluid to an elastic 
disordered solid, the gel. The viscosity coefficient grows 
with a power law behaviour, as function of the relative 
difference from the critical polymer concentration, char- 
acterized by a critical exponent k. The onset of the elastic 
response in the system, as function of the same control 
parameter, displays a power law increasing of the elas- 
tic modulus described by a critical exponent / [1]. This 
corresponds to the constitution inside the sol of a macro- 
scopic polymeric structure [2], that characterizes the gel 
phase. As implicitly suggested in the work of Flory and 
Stockmayer [2], the percolation model is considered as 
the basic model for the chemical gelation transition and 
the macromolecular stress-bearing structure in these sys- 
tems is a percolating network [3-5] . 
Moreover slow dynamics characterize the gelling solution: 
the relaxation functions in the experiments display a long 

time stretched exponential decay ~ e ^^o' as the gela- 
tion threshold is approached. The relaxation process be- 
comes critically slow at the gel point, where the onset of 
a power law decay is typically observed [6]. 
Gelation phenomena are also observed in colloidal sys- 
tems, that are suspensions of mesoscopic particles in- 
teracting via short range attraction: due to aggregation 
phenomena at low density (colloidal gelation) these sys- 
tems display gel states with a power law behaviour of the 
viscosity coefficient and of the elastic modulus [7] , as in 
chemical gelation. Yet the aggregation process gives rise 
to cluster-cluster aggregation producing a spanning clus- 
ter with a fractal dimensionality smaller than the ran- 
dom percolation case [8-11]. On the other hand with 
a weaker attraction at higher density a gelation charac- 



terized by a glass-like kinetic arrest [12,13] may be ob- 
served. The relaxation patterns closely recall the ones ob- 
served in glassy systems and are well fitted by the mode- 
coupling theory [14] predictions for supercooled liquids 
approaching the glass transition [13]. On the theoreti- 
cal side the application of the mode-coupling theory to 
systems with short range attractive interaction [15-17] 
(attractive glasses) has been recently considered and the 
connection with the colloidal glass transition has been 
proposed. 

The question, that we want to investigate in this paper, 
is whether and to what extent colloidal gelation, colloidal 
glass transition and chemical gelation are related and if 
a unifying description is possible. 

Let us first consider the case of the chemical gelation. By 
means of MonteCarlo numerical simulations on the cubic 
lattice we study a solution of tetrafunctional monomers. 
Each monomer occupies a lattice elementary cell and, to 
take into account the excluded volume interaction, two 
occupied cells cannot have common sites. At t = we 
fix the fraction (j> of present monomers respect to the 
maximum number allowed on the lattice, and randomly 
quench bonds between them. This actually corresponds 
to the typical chemical gelation process that can be ob- 
tained by irradiating the monomeric solution. The four 
possible bonds per monomers, randomly selected, are 
formed with probability pb along lattice directions be- 
tween monomers that are nearest neighbours and next 
nearest neighbours. Once formed, the bonds are perma- 
nent. The monomers diffuse on the lattice via random 
local movements and the bond length may vary but not 
be larger than Iq according to bond-fluctuation dynamics 
[20]. The value of Iq is determined by the self-avoiding 
walk condition and on the cubic lattice is = VTO in 
lattice spacing units. We let the monomers diffuse to 
reach the stationary state and then study the system for 
different values of the monomers concentration. We have 
considered pb = 1, for which the system presents a per- 
colation transition at 0^ = 0.718 ± 0.005. We have in- 
troduced this lattice model [18,19] to study the viscosity 
by means of the relaxation time, that diverges as func- 
tion of (j) sA, the percolation threshold with a power law 
behaviour [18]. The elastic response in the gel phase has 
been studied by means of the fluctuations in the free en- 
ergy and grows with a power law behaviour as well as 
function of {(j) — (j)c) [19]. In this letter instead we present 
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a completely new study of the dynamics in gelation phe- 
momena. 

To investigate the nature of the dynamic transition at the 

chemical gelation here we study the equilibrium density 
fluctuation autocorrelation functions fq-{t) given by 

<\p,<t')\>' 

where Pq{t) = X^ili e~*'^'''»(*\ ri(i) is the position of the 
i — th monomer at time t. N is the number of monomers 
and the average (...) is performed over the time t' . Due to 
the periodic boundary conditions the values of the wave 
vector q on the cubic lattice are q — ^{n,x,ny,nz) with 

= 1...L/2 integer values. 
We also study the mean square displacement of the par- 
tides (r^it)) = i Eliim + t') - n{t')?)- 
In Fig.l we present these; time autocorrelation functions 
as function of the time calculated on a cubic lattice of 
size L = 16. The data have been averaged over ~ 10 up 
to 10^ time intervals and over 20 different initial config- 
urations of the sample. As the monomer concentration 
(j) approaches the percolation threshold (/>c, fq{t) displays 
a long time decay well fitted by a stretched exponential 
law - e-f*/^)" with a /3 - 0.30±0.05. At the percolation 
threshold the onset of a power law decay is observed as 
it is shown by the double logarithmic plot of Fig.l with 
an exponent c [6]. As the monomer concentration is in- 
creased above the percolation threshold in the gel phase, 
the long time power law decay of the relaxation functions 
can be fitted with a decreasing exponent c, varying from 
c ^ 1. at (/)c to c ^ 0.2 well above 0c, where a nearly 
logarithmic decay appears. This behaviour well agrees 
with the one observed in gelling systems investigated in 
the experiments of refs. [6]. It is interesting to notice 
that this kind of decay with a stretched exponential and 
a power law reminds the relaxation behaviour found in 
spin-glasses [21]. 

The mean square displacement of the particles {r^{t)) 
presents a long-time diffusive behaviour and the diffu- 
sion coefficient decreases but remains finite also above 
<j)c. However the diffusion coefficient of clusters of size 
comparable with the connectedness length goes to zero 
at (pc with the same exponent as the relaxation time [18]. 
In colloids the aggregation is due to a short range at- 
traction and in general the monomers are not perma- 
nently bonded. To take into account this feature we 
introduce a novel ingredient in the previous model by 
considering a finite bond lifetime and study the ef- 
fect on the dynamics. It is worth noticing that the fea- 
tures of this model with finite Tf, can be realized in a 
microscopic model: a solution of monomers interacting 
via an attraction of strength —E and excluded volume 
repulsion. Due to monomers diffusion the aggregation 
process eventually takes place. The finite bond lifetime 
Tfc corresponds to an attractive interaction of strength 



—E that does not produce permanent bonding between 
monomers, and T(, ~ e^l^^ . A more detailed discussion 
will be available in a longer paper following up [22]. 
Due to the finite t;,, in the simulations during the 
monomer diffusion the bonds between monomers are bro- 
ken with a frequency 1/t{,. Between monomers separated 
by a distance less than ?o a bond is formed with a fre- 
quency /(,. For each value of t\, we fix /f, so that the 
fraction of present bonds is always the same [22] . 
After the system has reached the equilibrium we have 
calculated /g-(t) as defined in eq.(l) and the relaxation 
time T as /^(t) 0.1. 

Let us first analyze the behaviour of r, plotted in Fig. 2 
as function of the monomer concentration </> for different 
Tfe. For comparison we have also shown the behaviour of 
T in the case of permanent bonds, which displays a power 
law divergence at the percolation threshold 0c- We no- 
tice that for finite bond lifetime the relaxation time 
increases following the permanent bonds case (chemical 
gelation), up to some value 0* and then deviates from 
it. The longer the bond lifetime the higher 0* is. In 
the high monomer concentration region, well above the 
percolation threshold, the relaxation time in the finite 
bond lifetime case again displays a steep increase and a 
power law divergence at some higher value. This trun- 
cated critical behaviour followed by a glassy-like transi- 
tion has been actually detected in some colloidal systems 
in the viscosity behaviour [23,24]. 

These results can be explained by considering that only 
clusters whose diffusion relaxation time is smaller than 
Tf, will behave as in the case of permanent bonds. Larger 
clusters will not persist and their full size will not be rel- 
evant in the dynamics: the finite bond lifetime induces 
an effective cluster size distribution with a cut-off, which 
keeps the macroscopic viscosity finite [25]. As the concen- 
tration increases the final growth of the relaxation time 
is due to the crowding of the particles. 
In Fig. 3 /^(t) is plotted as function of time for a fixed 
value of Tb for increasing values of the monomer concen- 
tration (0 varies from below to well above 0c = 0.718). 
For small concentrations the autocorrelation function 
/q-(t) is well fitted by a stretched exponential decay, while 
for high monomer concentrations it c^xhibits a two-step 
decay, that closely resembles the one observed in su- 
percooled liquids. We fit these curves using the mode- 
coupling /3-correlator [14], corresponding to a short time 

power law ~ / -I- ^ \o^% time von Schweidler 

law ~ / — ^K^^ ' S^'^i^S exponents a ~ 0.33±0.01 
and h ~ 0.65±0.01 (the fuff lines in Fig.3). At long 
times the different curves obtained for different col- 
lapse into a unique master curve by opportunely rescal- 
ing the time via a factor t(0). The master curve is 
well fitted by a stretched exponential decay with /3 
0.50±0.06. The characteristic time r(0) diverges at a 
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value (j)g 0.96±0.05 with the exponent 7 ^ 2.3±0.1. 
This value well agrees with the mode-coupling prediction 
7 = l/2a+l/26 [14]. 

In the inset the data arc compared with the case of per- 
manent bonds for two different values of 0. We observe 
that for short time (of the order of Tft) the autocorrela- 
tion functions coincide with the permanent bond case. 
A similar pattern, characterized by a plateau, is found 
in the mean square displacement {f^{t)) (Fig. 4), with the 
diffusion coefficient going to zero as approaches 4)g . 
The inset also indicates that for short time the mean 
square displacement coincides with that calculated for 
the permanent bond case. 

These results show that when bonds are permanent 
(chemical gelation) the divergence of the relaxation time 
is due to the formation of a macroscopic critical cluster 
and the autocorrelation function exhibits a one step de- 
cay related to the relaxation of such cluster. In the case 
of finite Tf, there is an effective cluster size distribution, 
as discussed above, and the autocorrelation function ex- 
hibits a two step relaxation which is well fitted by the 
mode coupling theory. This can be explained if we con- 
sider the system made of effective clusters playing the role 
of single molecules in an ordinary supercooled liquid or in 
a colloidal hard sphere system. As the monomer concen- 
tration increases, the relaxation time will diverge due to 
the jamming of the effective clusters. The first decay is 
due to the rattling and relaxation of the " effective" clus- 
ter in a cage made of surrounding clusters. The second 
relaxation is due, as usual, to the opening of the cage re- 
sulting in a structural relaxation. This picture supports 
the jamming of clusters that has been suggested on the 
basis of experimental observations on colloidal gelation 
in ref. [12]. For high temperature t;, 0,the clusters 
reduces to single monomers with a crossover to the usual 
phenomenology of the hard spheres glass transition. 
When Tb is large enough (strong attraction) the cluster 
effect will dominate and the slow dynamics will exhibit 
features more closely related to chemical gelation (Fig.l). 
The only difference is that in the limit T5 ^ 00 we expect 
that the spanning cluster will have the structure of the 
cluster-cluster irreversible aggregation model instead of 
random percolation. 

In conclusion, these results suggest a unifying approach 
for chemical gelation, colloidal gelation and colloidal glass 
transition. In chemical gelation and colloidal gelation the 
cluster formation should produce the slow dynamics, that 
is expected to be of the same type of Fig. 1. In colloidal 
systems for weak attraction and high concentration the 
system crosses over from colloidal gelation to colloidal 
glass due to the jamming of effective clusters. 
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FIG. 2. The average relaxation time as function of the 
density; from left to right: the data for the permanent 

bonds case diverge at the percolation threshold with a 
power law (the full line); the other data refer to finite 
Ti, = 3000, 1000, 400, lOOMCstep/parttde decreasing from 
left to right (the dotted lines are a guide to the eye) 




FIG. 1. Double logarithmic plot of the autocorrelation 
functions /q-(t) as function of the time for q ~ 1.36 and 
= 0.6, 0.718, 0.8, 0.85. For (j> < (pc the long time decay is 
well fitted by a function (full Une) ~ e"^*/'^)'' with ~ 0.3. 
At the percolation threshold and in the gel phase in the long 
time decay the data are well fitted by a function ~ (1 + -^)"'^. 
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FIG. 3. The fq{t) as function of the time for 
q ~ 1.36 calculated on a cubic lattice of size L = 16 
for Tb = IQQOMC step /particle. (from left to right 
(j) = 0.6,0.718,0.8,0.85,0.87,0.9). The full lines correspond 
to the fit with the mode-coupling beta-correlator. In the in- 
set, the data for </> = 0.718 and = 0.8 (from left to right) are 
plotted together with the ones obtained at the same densities 
in the case of permanent bonds (dashed line). 
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FIG. 4. The mean-square displacement {f^(t)) of the 
particles as function of the time in a double logarithmic 
plot for Tb = WOOMC step /particle: from top to bottom 
(f) = 0.8,0.82,0.85,0.9, approaching (pg. In the inset (r^(t)> 
for the same 4> for rt, = WOOMC S /particle and in the 
case of permanent bonds (dashed line) (from top to bottom, 
= 0.718,0.8) 
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